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Abstract Achieving sub-arcsecond co-registration across varying time-lines of 
multi-wavelength and instrument images is not trivial, and requires accurate 
characterization of instrument pointing jitter. In this work we have investigated 
internal pointing errors, on daily and yearly time-scales, occurring across the 
Solar Dynamics Observatory's (SDO) Atmospheric Imaging Assembly (AIA) 
and Helioseismic Magnetic Imager (HMI). Using cross-correlation techniques 
on AIA 1700 A passband and HMI line-of-sight (LOS) magnetograms, from three 
years of observational image pairs at approximately three day intervals, internal 
pointing errors are quantified. Pointing variations of ±0.26" (jitter limited) and 
± 0.50" in the solar East-West (x) and North-South (y) directions, respectively, 
are measured. AIA observations of the Venus June 2012 transit are used to 
measure existing coalignment offsets in all passbands. We find AIA passband 
pointing variations are (A Xqo) = 1.10" ±1.41" and (A Yqo) = 1.25" ± 1.24", 
when aligned to HMI’s nominal image center, referred to herein as the CutOut 
technique (CO). Minimal long-term pointing variations found between limb 
and correlation derived pointings provide evidence that image center positions 
provided by the instrument teams achieve single pixel accuracy on time-scales 
below their characterization. However, daily AIA passband pointing variations 
of <1.18" indicate autonomous sub-arcsecond co-registration is not yet fully 
achievable. 


1. Introduction 


Addressing the coronal heating problem is not trivial, as it requires the use 
of high resolution multi-wavelength and temporal images with sub-arcsecond 
co-registration. Currently, in the solar physics field the Solar Dynamics Ob¬ 


servatory's (SDO) Atmospheric Imaging Assembly (AIA; Lemen et al. 2012) 
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and Heliographic Magnetic Imager (HMI; Schou et al. 2012) instruments are 
providing unprecedented amounts, approximately one image a second, of high 
resolution data (<0.6"). AIA takes full-disk images of the Sun in the following 
nine passbands: 94 A (logT ss 6.8), 131 A (logT « 5.8), 171 A (logT ss 5.9), 193 A 
(logT«6.2, 7.2), 211A (log T« 6.3), 304A (logT«4.8), 335A (log T« 6.4), 
1600 A (logTss 5.0), 1700 A (logT~ 3.7), and 4500 A (logTss 3.7; observed typ¬ 
ically every ss 30 min), while HMI takes full-disk images of the Sun’s line-of-sight 
(LOS) magnetic field (~45 s). 

Some recent studies have highlighted that commonly used multi-passband 
alignments for AIA and HMI, considered to yield sub-pixel (< 0.6") accuracy, de- 


liver a less than optimized status (s 

M";I 

3el Zanna, O’Dwyer, and Mason 

2011 

Brooks, Warren, and Ugarte-Urra 

201 

2). These alignment techniques, 

both 


of which utilize aia.prep.pro, i.e., standard Solar Software (SSW), and result 
in level-1.5 data, are summarized as follows: aligning multi-passband charge- 
couple device (CCD) centers as Sun center (SC) or using a single passband’s 
documented alignment information as a fiducial reference position for all other 
passbands. Thermal jitter motion, which places a systematic limit on achievable 
co-alignment accuracy, reported to affect each AIA telescope is K3 0.3" (Lemen 


et al, 2012 

Aschwanden et al., 2013 

1, while that of HMI remains un-reported 

(as of time of writing; 

Schou et al. 

2012 

)• 


Shine et al. (2011) showed AIA passband image pointings are not temporally 


stable, suggested as a result of thermal flexing, for the four AIA telescopes and 
indicated monitoring and characterization of each telescope’s nominal pointing 
(i.e., offset relative to solar center) is required. Currently, per AIA passband, 
image pointings are defined from an automated limb finder routine that performs 
an interactive Hough transform on candidate limb pixels to converge on the solar 
center (x, y) and define it’s radial distance (Shine et al, 2011). A daily average 
of these image center positions, performed once a week, are used to define the 


master pointings that are propagated in documented header information (Shine 


et al., 2011). However, image positions change significantly on weekly time scales, 


and daily variations (>0.6") are not characterized by master pointings (Shine 


et al., 20111. AIA image plate scales and angles are considered significantly 


more stable than that of the aforementioned offsets, with Shine et al. (2011) 
indicating no evidence of scale change and angle variations of 0.001°. HMI’s 
master pointings are defined from similar techniques and time scales as those 
used by AIA. 

Using cross-correlation techniques on AIA and HMI observational image pairs, 
performed on sub-weekly time-scales, we measure instrument jitter in the master 
pointing data by directly comparing limb and correlated image center positions. 
The June 2012 Venus transit is used to measure pointing variations of each of 
AIA’s passbands. Throughout the remainder of this paper the term SC tech¬ 
nique refers to multi-passband image CCD centers defined as Sun center, while 
CutOut is that where HMI’s master pointing is as the fiducial reference position 
for all AIA passbands. It is noted, reported co-registration discrepancies are 
consistent across both co-alignment techniques summarized above, but result in 
a constant solar center offset « 5" (Figure [l]) between the two (Boerner 2013, 
private communication). 
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Table 1. Observational date range, image time differ¬ 
ence ((At); in seconds), typical observed x and y image 
center offsets ((i?J bs ) and (R° hs ), respectively, in arc- 
sec) for co-aligned 1700 A far-ultraviolet (FUV) images 
and LOS magnetograms. 


Obs. Date Range 

(At) 

(R% bs ) 

(R° hs ) 

May 2010 - March 2013 

5.71 

0.26 

0.50 


In that respect, the remainder of this paper is outlined as follows: in Sec- 
tion [2] we provide a long-term, sub-weekly, characterization of AIA and HMI 
instrument jitter and compare image pointings resultant from limb fitting and 
cross-correlation techniques; Section [3] quantifies AIA passband pointing varia¬ 
tions via observational data from the June 2012 Venus transit; a discussion of our 
findings is provided in Section [4j and concluding remarks are given in Section [5] 


2. AIA and HMI Sub-Weekly Pointing Errors 


Measurements of AIA and HMI’s instrument pointing variations are performed 
at approximately three day intervals spanning May 2010-March 2013 from 
cross-correlated image pairs of the 1700 A passband (which directly observes the 
solar photosphere) and LOS magnetograms. Cross-correlation, widely considered 
to yield sufficient accuracy (< 1 pixel; e.g., Brooks, Warren, and Ugarte-Urra 


2012J Del Zanna, O’Dwyer, and Mason 2011) is done using visually bright 

point-like features or magnetic neutral lines (Figure [l]). 

Image pairs were pre-processed to level-1.5 using aia_prep.pro (which included 
roll angle corrections), CutOut aligned, and magnetograms interpolated to the 
effective resolution of FUV i ma ges (~ 0.6"). Solar rotation effects, «0.17 arcsec 
min -1 (Sarro and Berihuete, 2011) close to disk center, were minimized by 


using observational time differences below AIA’s thermal jitter motion (i.e., 
< 2 min), and the standard SSW routine drot_map.pro. FUV and magnetogram 
images were then mapped to solar coordinates by calibrating each image’s axes 
with respect to actual Stonyhurst heliographic coordinates, i.e, coordinate (0,0) 
corresponds to solar center. Disk regions of ± 300" from solar center are only 
considered to minimize projection effects. Per observational image pair, two 
measurements of the x and y co-alignment offsets, relative to HMI’s master 
pointings (l?° bs and i?° bs , respectively), were made on regions comprised of 
bright network and small compact features {e.g., similar to those observed in 
Figure [lj. The resultant l?° bs and R° hs offsets were then smoothed per day with 
errors propagated from the differences in daily measurements. 

Figure [ 2 ] displays AIA and HMI’s master pointings {i.e., XC/YC mp derived 
from limb fitting techniques) compared to the pointings obtained from cross- 
correlating image pairs (i.e., XC/YC 0 b s ) at sub-weekly time-scales. A summary 
of our co-alignment offsets are provided in Table [l] with typical values of ss 0.3" 


SOLA: vtstudy.tex; 22 January 2015; 1:26; p. 3 
















Orange et al. 


15 May 2010 06:05:59 UT 




\ ^ 


SC 


15 May 2010 06:05:59 UT 

s' 


X 


CutOut 

1 15 May 2010 06:05:59 UT 


Figure 1 . The top and middle panels (pre-processed under SC and CutOut techniques, respec¬ 
tively) are AIA 1700 A intensity images (DN pix —1 s —1 ), shown on inverted color scale, with 
contours showing regions of positive magnetic flux (blue solid lines) derived from HMI’s LOS 
magnetograms and plotted at levels of zb 25-100 G. Bottom panel is AIA 1700 A same as that 
in the top but contours are from the middle panel image, which shows the difference between 
the pre-processing techniques is a solar center offset (« h" in both the x and y directions). 


and 0.5 7/ (in solar x and y directions, respectively). These results are consistent 
with reported mis-alignments (< 1"; Del Zanna, O’Dwyer, and Mason 2011 


Brooks, Warren, and Ugarte-Urra 2012) as well as those expected from daily 
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Figure 2. Top, from left to right, are plots of AIA and HMI master pointings, respectively, 
(XC mp ; asterisks) provided by the SDO teams and observed image centers from cross-corre¬ 
lation techniques ( AC 0 bsl red pluses) in the solar x direction versus time (May 2010-March 
2013), respectively. Bottom, similar to top panels for pointings in the solar y direction. 


variations (< 1.2"; Shine et al. 2011). Moreover, i?° bs is within the limits of 


instrumental jitter, and as observed in Figure [2] is indicative of random scat¬ 
ter. Solar y pointing offsets varied minimally (ss ± 0.02") over the three year 
span, and indicate their alignment agreement could be autonomously improved 
(Figure [2]). 


3. Using the Venus Transit for Investigation of Internal AIA 
Co-Registration Errors 


Investigation of daily AIA pointing variations uses observational data of the 5-6 
June 2012 Venus transit (23:00:00 UT-04:00:00 UT), while Venus’ silhouette 
was fully visible on the solar disk (Figure j3j. Venus’ transit provides the ideal 
fiducial reference for measuring such errors since its silhouette must be seen 
at identical locations in each passband when the transit speed across the solar 


surface is accounted for (Kamio et al., 2010 Shimizu et al., 2007). 


Daily pointing variations from the various AIA passbands are determined 
from a direct comparison of the heliographic coordinates of Venus’ center to those 
predicted by deriving the Venus’ velocity using running difference images (1700 A 
passband; see Figure [4]). The details of our observational data and method as 
well as results are described in the next two Sections 3.1 and 3.2. 
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Figure 3. AIA 193 A passband observed on 6 June 2012 at 00:10:01 UT. The image is 
displayed on a log 1(l scale and the boxed blue region surrounds the disk of Venus. 


3.1. Venus Transit Analysis 

The Venus transit data were pre-processed to obtain both SC and Cut-Out multi¬ 
image alignments, as well as those aligned to their respective master pointings 
(MP). For each pre-processing technique, a ten minute temporal cadence was 
used which resulted in « 30 images of each passband with exception of 4500 A 
images (5 observations). Custom written software, discussed in depth below, was 
then used to measure the heliographic coordinates of Venus’ disk centroid from 
all images. 

The software measured individually the vertical, North - South (N-S), and 
horizontal, East - West (E-W), disk centroid positions vi a calculating and 
comparing observed disk diameters to the planet’s known angular size (58.3"; 
Odenwald 2012). Calculation of the disk diameter required user selection of loci 


regions corresponding to a single orientation’s limb position ( e.g ., leading and 
trailing edges for E-W, see Figure [5]). Per loci region, limb positions along the 
vertical or horizontal for E-W and N-S diameters, respectively, were determined 
by measuring where the steepest flux gradients occurred (Figure[5]). The position 
of steepest flux gradient occurring furthest away from disk center was stored and 
visually checked. Uncertainties in each of the aforementioned measurements were 
derived from the plate scale of a single pixel per passband (defined in the header 
of the FITS data file) and residuals from comparing measured angular sizes to 
the known value. 
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Figure 4. Sequence of AIA 1700 A images, observed on 5-6 June 2012 from 23:00:07 UT 
- 03:52:07 UT (left to right, respectively) with pluses denoting Venus’ centroid, defined by 
the technique discussed in the text (Section |3.1| ). Arrows are shown to represent the running 
difference technique applied to the data to measure the planet’s transit velocity in the x and 
y directions and provide fiducial references as a function of time from which co-alignment 
discrepancies between all other AIA passbands are measured. 


; P- 7 
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Figure 5. The top panel is an intensity image centered on the silhouette of Venus’s disk from 
AIA’s 1700 A passband observed on 5 June 2012 23:20:07 UT, with the center of the crosshairs 
representing the measured center of Venus’ disk. The four bottom panels are flux versus pixel 
position for the North, East, South, and West positions noted in the top panel from top left to 
bottom left in clockwise direction, respectively. Note, vertical dotted blue lines denote measured 
limb positions, as outlined in the text, of (-385.4", 673.5^), (-385.4 /; , 673.5 ,/ ), (-355.5 /; , 644.5 // ), 
and (-413.5 7/ , 644.5 ;/ ) for the North, East, South, and West positions, respectively, presented 
here in the form of solar x and y directions, respectively. 


Figure [5] provides an example of this entire process. Its left panel displays 
the resultant diameters (solid lines, plotted with a width of «1" to be easily 
observable) and limb positions (marked by their respective direction and a solid 
line). The right panel displays flux versus position for the N-S and E-W disk 
edges corresponding to the aforementioned marked locations on the intensity 
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Figure 6. Intensity images centered on the silhouette of Venus’ disk from AIA’s 1600 A, 304 A, 
94 A, and 131 A passbands from top right to bottom right in clockwise direction, respectively, 
observed on 5-6 June 2012 showing variations in the contrast of Venus’ disk as observed by 
different passbands. On each image the center of the crosshairs represents Venus’ measured 
centers determined from the techniques discussed in the text, respectively. 


image of the right panel. Figure [6] has also been provided to illustrate the 
variations in Venus’ contrast against the solar disk for a sample of passbands. 
It is noted, features which show up on Venus’ disk in 131 A and 94 A images 
of Figure [6] are directly indicative of the stray light which affects all AIA EUV 
passbands (Poduval et al ., 2013). 

Using running differences of Venus’ disk center position, its transit velocity 


1-4 


,-i 


s ' ana v y = — y.^x IU _ 5 ± 6.3x 10 -4 arcsec s -1 in the horizontal and vertical 
directions, respectively (Figure [4]). Observational disk centers and the transit 
speed were then used to generate a set of predicted disk positions at one second 
intervals over the observational time frame used. Co-alignment offsets (AX, 
AY), relative to 1700A images as function of transit time, were measured for 
co-aligned multi-passband data. 
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Table 2. Typical co-alignment offsets, in arcsec, 
(averaged over seven AIA passbands, and mea¬ 
sured relative to the 1700 A image positions) and 
their resultant uncertainties, defined as the stan¬ 
dard error on the mean (SEM), for pre-processing 
images under SC, MP, and CutOut assumptions. 


Pre-Process 

(X) 

a SEM 

(Y) 

Y 

a SEM 

sc 

1.47 

2.86 

1.83 

2.29 

MP 

1.29 

1.60 

1.37 

1.66 

Cut-Out 

0.87 

1.14 

1.18 

1.46 


3.2. Results 


Figures [7] and [s] provide the x and y co-alignment offsets, respectively, (plotted at 
ss 30 min intervals for clarity) of AIA passbands imaged by telescopes 1, 2, and 4 
( |Lemen et al., 2012| . Figure [9] displays both the x and y co-alignment offsets of 
AIA passbands imaged by telescope 3 ( |Lemen et al., 2012 1. As observed in these 
figures, significant variance exists, not only between disk center positions but 
as a function of time, too. Most typically, the x and y pointing variations are 
indicative of random scatter, but a number of passbands do exhibit a quasi- 
periodic nature in the x direction. Using the 304 A passband the period of 
this variation was estimated to be sa 5 h. However, we note that it is only a 
crude approximation since variations are found in other passbands (~±30 min, 
Figures pf] and |9| while our observation sequence did not cover a full period. 

Table|2] presents the co-alignment offsets, averaged over seven passbands, and 
uncertainties, defined as the standard error on the mean (SEM), for each pre¬ 
processing technique. It is noted, these results are not weighted by measurements 
derived from 94 A and 4500 A images f or the f ollowing reasons: the low count- 
rates of the 94 A passband (Boerner et al., 2012), and the significantly longer 
temporal cadence of the 4500 A passband. The results of Table [ 2 ] indicate that 
image pointings vary by > 2 pixel, independent of solar direction, on daily time- 
scales. We point out, no technique results in sub-pixel (<0.6") alignment, as 


expected (Shine et al., 2011), while resultant uncertainties indicate minimal 


variations of ss 1.3" occur over a few hours. 

Co-alignment offsets per passband were re-measured by letting individual 
passbands serve as the fiducial reference (i.e., using running image differences 
to derive and predict the planet’s disk center positions), again barring 94 A and 
4500 A images. The results of this analysis are shown in Table [ 3 ] As observed, 
the typical alignment accuracy, as well as its respective uncertainty, is consistent 
with our previous reports for the 1700 A passband. Therefore, using our resultant 
co-alignment offsets, the typical AIA passband alignment accuracy, derived from 
eight passbands, is 


(A X CO ) = 1"10±1"41 and, 
(A Y co )= 1"25± l."24, 
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Figure 7. Solar x co-alignment offsets (AX), relative to 1700 A image positions, versus time 
(minutes measured relative to 5 June 2012 23:00:00 UT) shown as a function of AIA telescope 
number (1, 2, and 4 only; Lemen et al., 2012) and smoothed over £3 30 min for visual clarity. 
Pre-processing techniques of MP, SC and CutOut are represented by red pluses, blue x-signs, 
and gold triangles on each plot, respectively. 
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Figure 8. Solar y co-alignment offsets (AY), relative to 1700 A image positions, versus time 
(minutes measured relative to 5 June 2012 23:00:00 UT) shown as a function of AIA telescope 
number (1, 2, and 4 only; Lemen et al., 2012) and smoothed over p^ 30 min for visual clarity. 
Pre-processing techniques of MP, SC and CutOut are represented by red pluses, blue x-signs, 
and gold triangles on each plot, respectively. 
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Figure 9. Solar x and y co-alignment offsets (AA' and AY in left and right columns, respec¬ 
tively), relative to 1700 A image positions, versus time (minutes measured relative to 5 June 
2012 23:00:00 UT) shown for AIA telescope number 3 (Lemen et al., 2012) and smoothed over 
Ri 30 min for visual clarity. Pre-processing techniques of MP, SC and CutOut are represented 
by red pluses, blue x-signs, and gold triangles on each plot, respectively. 


in the x and y directions, respectively, independent of the reference band. 


4. Discussion 

Using an extended three year study, a characterization of AIA (1700 A) and 
HMI’s (LOS magnetograms) pointing variations, at sub-weekly time-scales, was 
performed. Mis-alignments of the two instruments were consistent with previous 
reports indicating no better than single pixel accuracy exists in the master 


pointing data 1 

Del Zanna, O’Dwyer, 

and Mason, 2011 

Brooks, Warren, and 

Ugarte-Urra, 2012 

Shine et al., 2011 

). Pointing variations of the various AIA 


passbands were determined to be < LIS" on a daily scale using the recent June 
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Table 3. Typical co-alignment offsets, in arcsec, 
(averaged over 7 AIA passbands) measured rela¬ 
tive to the waveband listed in the far left column’s 
image positions, for pre-processing under CutOut 
assumptions only. 


Band A 

(Xco) 

a SEM 

(Yco) 

Y 

a SEM 

131 

0.95 

1.19 

1.24 

1.56 

171 

0.98 

1.27 

1.14 

1.39 

193 

0.96 

1.21 

1.19 

1.47 

211 

1.06 

1.31 

1.27 

1.63 

304 

1.31 

1.62 

1.41 

1.72 

335 

0.97 

1.22 

1.56 

1.89 

1600 

1.73 

2.32 

1.08 

1.32 


2012 Venus transit, and as such were consistent with the instrument team’s 
expectations (<1.2" 


Shine et al. 2011) 


Two notable artifacts were identified in the long-term pointing data (Fig¬ 
ure]^. The first was a significant shift in AIA’s pointing between December 
2011 and January 2012, while the second was annual sinusoidal fluctuations in 
both HMI and AIA’s master pointings (Figure [2|. The final artifact was daily 
variations found in the master pointings of AIA passbands (Figures [7j |s| and [9]). 
Below we provide discussions on the causes of these two artifacts as well as how 
they relate to the information obtained in this article. 

The AIA 1700 A pointing shift (Figure [ 2 ]) is due to adjustments made to 
telescope 3’s thermal control parameters to increase the instrument stability 
(Boerner 2013, private communication). The heater adjustment occurred on 18 
January 2012 and was also performed on telescopes 1 and 2. These types of 
internal updates have been performed at various times in the mission, and in 
each case updates are made to master pointings to compensate for the resultant 
shifts. 

The seasonal drift in telescope pointings is likely caused by thermal flexing of 
the instrument. This flexing is considered to occur between the science telescope 
(i.e., makes passband images), and the guide telescope which drives the image 
stabilization system (Boerner 2013, private communication). Changes in energy 
flux due to variations of SDO’s distance to the Sun during its yearly orbit, with 
peak changes correlating with Earth’s aphelion and perihelion, are suggested 
as the origin of the flexing itself. Furthermore, SDO’s geosynchronous orbit 
likely generates these same effects on a daily basis, thereby causing similar daily 
variations in the master pointings. Temperature variations are believed to be 
the cause of daily pointing variations («0.6" 1.2"; Shine et al. 2011), to which 
our results are consistent (Table [3|. 

AIA and HMI long-term pointing variations in the solar x direction, «0.26", 
indicates the master pointing accuracy is jitter limited at time-scales below their 
characterization. This conclusion is further supported by the distribution ( i.e ., 
random scatter) of our measured x pointings to those of the instrument teams 
(Figure [ 2 ]). It is hypothesized that the consistent nature of the y offsets, between 
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limb and correlation derived pointings, may be an indication of expected satellite 
flexing or even a payload shift. Finally, based on minimal variations of correlated 
offsets between AIA and HMI we hypothesize the thermal jitter affecting HMI 
is similar in scale to AIA’s. 


5. Conclusions 


Our study has provided the first characterization of internal pointing errors 
occurring across SDO’s AIA and HMI instruments using cross-correlation tech¬ 
niques at time scales below the team’s master pointing data. Long-term cor¬ 
relation offsets provided an estimate of pointing jitter affecting the AIA and 
HMI instruments. Our solar x pointing variations («0.3 ,/ ), between limb and 
correlation techniques, are within the limits of thermal jitter in the AIA in¬ 


strument (Lemen et al. , 2012 Aschwanden et al. , 2013), while those of the 


solar y direction («0.5 ,/ ) are consistent with expected variations (Shine et al., 


2011). Thermal jitter, as previously highlighted, occurs from thermal bending 


between the science and guide telescopes as a result of orbit induced temperature 
variations that occur on yearly and daily time scales. Minimal variations between 
limb and correlation derived pointings indicate that master pointing data achieve 
single pixel accuracy on sub-weekly time-scales. Long-term instrument pointing 
variations between HMI and AIA also suggest that thermal jitter affecting HMI 
is similar in scale to that of AIA. 

From the June 2012 Venus transit, daily pointing variations in all AIA pass- 
bands were presented. Our results are consistent with the instrument team’s 


expectations of daily pointing shifts of 1-2 pixels (Shine et al., 
typical co-alignment offsets were (AAVo) = 1-10" 


2011). The 
when 


and (AYco) = 1-25' 
aligned to HMI’s nominal image center. Finally, these co-alignment accuracies 
are consistent with previously reported mis-alignments of ss 1" in AIA passbands 


(Del Zanna, O’Dwyer, and Mason, 2011 Brooks, Warren, and Ugarte-Urra, 


2012), and expected variations in master pointing data (Shine et al., 2011). 
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